First- and second-order superfluid-Mott-insulator phase transitions of spin-1 bosons with coupled 

ground states in optical lattices 
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We investigate the superfluid-Mott-insulator quantum phase transition of spin- 1 bosons in an optical lattice 
created by pairs of counterpropagating linearly polarized laser beams, driving an Fg = 1 to _Fe = 1 internal 
atomic transition. The whole parameter space of the resulting two-component Bose-Hubbard model is studied. 
We find that the phase transition is not always second order as in the case of spinless bosons, but can be first 
order in certain regions of the parameter space. The calculations are done in the mean-field approximation by 
means of exact numerical diagonalization as well as within the framework of perturbaton theory. 
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PACS numbers: 03.75.Lm,03.75.Mn,71.35.Lk 

The superfluid-Mott-insulator quantum phase transition 
(SMQPT) of spinless bosons in periodic lattices is a second- 
order transion, which is characterized by a continuous varia- 
tion of the order parameter ijj from V' 7^ (superfluid phase) 
to -0 = (Mott-insulator phase) if the amplitude of the lattice 
potential increases |1, 2]. In our recent paper (Jl, we have 
shown that the SMQPT in a system of spin-1 bosons can be 
first order as well |4]. By means of numerical calculations 
within the framework of the mean-field theory, it was found 
that in the case of ^'^Na the SMQPT is second order if the 
number of atoms per lattice site n — 1,3, and it is first order 
for n = 2. In the case of ^''Rb, the SMQPT was found to be 
second order for n = 1, 2, 3. In the present work, we con- 
tinue the study of Ref. The main purpose is to investigate 
the whole parameter space of spin- 1 bosons and to find the re- 
gions where the SMQPT is first and second order for arbitrary 
n. 

We consider spin-1 neutral polarizable bosons, possessing 
three Zeeman-degenerate internal ground and excited elec- 
tronic states characterized by the magnetic quantum number 
m — Q, ±1 {Fg = i^e = 1) in a d-dimensional {d = 1, 2, 3) 
optical lattice. The lattice is assumed to be created by d pairs 
of counterpropagating linearly polarized laser waves running 
in d orthogonal directions and having different frequencies in 
different directions. The beams propagating along the 3-axis, 
which is chosen to be a quantization axis, are polarized along 
the 1- or 2-axis, and the beams propagating along the 1- and 
2-axis are polarized along the 2- and 1-axis, respectively. 

The running laser waves form left- and right- 
polarized standing waves with Rabi frequencies 
^uif,^) = ^loi^ cos{kLr^), which couple internal ground and 
excited states by V- and A- transitions. In order to avoid 
decoherence due to spontaneous emission, the detunings 
Ai must satisfy the conditions jA^j ^ 7, where 7 is the 
spontaneous emission rate. If the laser intensities and the de- 
tunings are chosen in such a manner that f^Q^/A^ = f^o/^' 
V = 1 , . . . , d, the laser potential acting on the atomic ground 
states is given by the matrix 
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riod tt/Zcl and at the same time couples the atomic ground 
states with m — ±1. In the case of red detuning A < 0, 
the system is described by the two-component Bose-Hubbard 
Hamiltonian 0] 
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where /i is a chemical potential and the operator fii — fiQi + 
fiM- The indices «, j label the lattice sites, and a^i is the Bose 
annihilation operator for the component cr = 0, A attached to 
the ith lattice site. The tunneling matrix element 
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where Wi{r) is the Wannier function for the lowest Bloch 
band of the potential 



yLo(r) = 2?i^^cos2(fcir,), 

is a rapidly varying function of the laser intensity ^9]. The 
variations of the atomic interaction parameters 



Us,a = gs,aJw,Hr)dr 



are much slower |6]. The quantities gs^a describe the repul- 
sive interaction of the condensate atoms and the spin-changing 
collisions. The parameter Us is positive, but Ua can be either 
positive or negative depending on the sign of ga- 
in order to investigate the SMQPT, we employ the mean- 
field approximation LL2J 
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which determines the isotropic lattice potential with the pe- 



where ij^a is the order parameter for Bose-Einstein conden- 
sation in the component a — 0, A, which can be considered 
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as a real quantity. In this approximation, the Bose-Hubbard 
Hamiltonian becomes local and every lattice site is described 
by the Hamiltonian 



(4) 



ij(0) = 2dJ {^l + ^l) + -fn{n - 1) 



+ Ua 

V = ~2dJ 



where the index 1(2) corresponds to Ua < (Ua > 0). In 
Eq.Q, we have omitted the site index i and introduced the 
isospin operator T with the components 

Ti = (^a]yaa + dl&A^ /2 , 

2^2 = i (d\dQ - aJoA^ /2 , (5) 

fa = (cilao - d^'^aA^ /2 , 

which has the property 
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If the tunneling V is negligible, the eigenstates of the 
Hamiltonian @ are determined by H^^'' and can be calculated 
analytically. They are eigenstates |n/2, A^) of the isospin op- 
erator with the corresponding eigenenergies given by 



E 



(0) 

n/2,M 



2dJ (V'o^ + ^2 ) ^ ^^(^ 



1) 



lin 



(6) 



where n is the number of atoms per lattice site and Ai ~ 
~n/2, . . . , n/2 is the isospin projection on the direction 1(2) 
in the isospin space for the case of Ua < 0{Ua > 0). Note that 
n/2 plays the role of the isospin quantum number If n is even, 
the states with Al = are not degenerate, while the others 
with M 7^ are doubly degenerate. If n is odd, the states 
with = do not exist and, therefore, all the eigenstates 
of iJ^") are doubly degenerate. If Ua < 0, the ground-state 
energy _Eg(V'o, V'a) — ^1^/2 ±ri/2' chemical potential 

/i varies in the interval 



1 < 



Us - \Ua\ 



< n 



n= 1,2,... 



In the case of Ua > 0, Egiijjo, Va) = -E'i/2,±i/2 ^^'^ 

Us{n-1) < I.K UsU-Ua , 



if n is odd. If Ua > and n is even, Eg{'>pQ, ip^) — E^^j^ ^ 
and 

Us{n -l)-Ua< H<Usn. 
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FIG. 1: Ground-state energy of the Hamiltonian ^4} for ^^Na 
(Ua/Us ^ 0.037 (J). p/C/« = 1.5, 2dJ/Us = 
0.125(1), 0.148(2), 0.157(3), 0.167(4). 



The minimum of Eg in all the cases is reached at -00 — V-'A — 
and the system has a vanishing compressibility; i.e., it is 
in the Mott-insulator phase. Since {n/2, M\f3\n/2, M) = 
0, if Ai is the isospin projection on the direction 1 
or 2, the mean occupation numbers of the componets 
{n/2,M\no\n/2,M) = {n/2,M\hA\n/2,M) ^ n/2. 

In order to investigate the transition from the Mott phase 
into the superfluid phase, one has to calculate the ground-state 
energy of the complete Hamiltonian 0. This can be done 
exactly by numerical calculations or approximately treating 
the tunneling as a perturbation. Let us consider the case 
Ua < first. Nonvanishing matrix elements of the operator V 
are given by 
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= -dJ^n ± 2M (7) 
X (?/;A±-(/'o)exp(=Fi7r/4), 



which implies that the ground-state energy Eg is indepen- 
dent of the sign of t/i^, £:g(-0A,V'o) = £^s(I^a|, iV'ol), and 
it is a symmetric function of and ■(/'a: Eg{ipo,4'A) = 
Eg{il)A,'4'o)- The minima of Eg are located on the lines 
i/^a — i'/'o f3'l. Therefore, for the calculation of the phase 
diagram of the system one can set ?/'a — ±V'o = '/'■ Employ- 
ing the perturbation theory we found that Eg has the structure 



Eg{il)) = ao + a2i/'^ + 04 V'^ + 



(8) 



where a2 can be either positive or negative, but a4 is always 
positive, which means that the SMQPT is second order. The 
boundary between the Mott and superfluid phases is deter- 
mined from the condition 02 = and given by 



J (1 - n + /t) (n - /t) / (1 + /i) 



(9) 



where J — 2dJ / {Us — \Ua\) and li — ^/ {Us — \Ua\), similar 
to the case of spinless bosons 1 1, 2]. This equation agrees per- 
fectly with the results obtained by numerical diagonalization. 
In the case of positive Ua, the situation is quite different. 
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FIG. 2: In the shaded regions of this diagram, 04 is negative and 
Eg{^) has two minima at certain values of J. In the remaining part 
04 is positive and Eg{tp) has only one minimum. 
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FIG. 3: Phase diagram for ^^Na (Ua/Us ^ 0.037 |7]). The regions 
of metastable superfluid phase coexisting with the stable Mott phase 
and that of metastable Mott phase coexisting with the stable super- 
fluid phase are denoted by MS and MM, respectively. 



Since the matrix elements of the operator V are given by 
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the ground-state energy Eg depends only on — i/jq + ipj^. 
Typical dependences Eg{ip) obtained by numerical diagonal- 
ization are shown in Fig. [2 In contrast to the case of negative 
Ua, there can be one or two minima. This leads to the fact 
that the QPT can be not only first order but also second or- 
der and the superfluid and Mott phases can coexist in certain 
parts of the phase diagram. Analytical calculations within the 
framework of the perturbation theory show that Eg has the 
form (|8j, but now not only 02 but also 04 can be either pos- 
itive or negative. Therefore, in order to be able to work out 
the complete phase diagram analytically, one has to calculate 
flg in Eq.(|8}. This is a tedious task and we have not done 
that. However, solving the equation a4{ii/Us, Ua/Us) ~ 
one can determine the regions in the plane [ji/Ug, Ua/Us), 
where Eg can have one and two minima; i.e., one can find the 
regions where the SMQPT is second and first order, respec- 
tively. These regions, which are shown in Fig.|2] can be also 



found by numerical diagonalization and the result appears to 
be the same. Our analytical as well as numerical calculations 
show that, for n = 1, Eg {ip) has only one minimum and the 
SMQPT is always second order If n > 2, Eg{ip) can have 
either one or two minima, depending on the parameters, and 
the SMQPT can be either second or first order. If, at a fixed n, 
Ua/Us is larger than some critical value, which is about 0.188 
for even n and grows from 0.012 (n — 3) to 0.015 (n — > 00) 
for odd n, the QPT is second order; otherwise, it is first or- 
der In the case of ^"^Na shown in Fig.Ol an interesting regime 
is achieved, when the QPT for odd n is second order, but for 
even n it is first order 

The phase diagram for ^'^Na presented in Fig. |3] has been 
obtained by numerical diagonalization. It consists of a series 
of (internal) lobes corresponding to the stable Mott phase and 
external regions corresponding to the stable superfluid phase. 
However, in the case of even n, the two regions are separated 
from one another by intermediate ones, where the stable and 
metastable superfluid and Mott phases coexist. The boundary 
separating the region of the stable superfluid phase from other 
ones can be determined from the condition 02 = 0. If n is 
odd, it is given by 
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n — U'a — n' /i' — n + 1 



,(11) 



where J' — 2dJ/Us, U'a — Ua/Us, and ^' — /i/C/s- For even 
n, the boundary is given by the equation 



J' 



{li'-n+l + U'a) jn - y.') 
11' -n + l + U'a + {l + U'a)n/2 



(12) 



These analytical expressions describe perfectly the corre- 
sponding numerical results. 

In conclusion, we have shown that in the laser configuration 
we consider, the SMQPT of spin-1 bosons with the ferromag- 
netic interactions (Ua < 0) is always second order and the 
boundary between the superfluid and Mott-insulator phases is 
given by Eq.Q. In the case of antiferromagnetic interactions, 
the SMQPT is first order in the shaded region of the plane 
(/i, Ua), shown in Fig.|2] and it is second order in the rest part 
of the plane. 
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